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We generalize a semiclassical theory and use the argument of angular momentum conservation 
to examine the ballistic transport in lightly-doped Weyl semimetals, taking into account various 
phase-space Berry curvatures. We predict universal transverse shifts of the wave-packet center in 
transmission and reflection, perpendicular to the direction in which the Fermi energy or velocities 
change adiabatically. The anomalous shifts are opposite for electrons with different chirality, and 
can be made imbalanced by breaking inversion symmetry. We discuss how to utilize local gates, 
strain effects, and circularly polarized lights to generate and probe such a chirality Hall effect. 


When a strong topological insulator [1, 2] undergoes 
a phase transition to a trivial band insulator in three 
dimensions, a gapless Dirac point [3-6] emerges at the 
critical point, if both time-reversal (T) and inversion 
(V) symmetries are present. Remarkably, when one of 
the two symmetries is broken, the critical point expands 
in the phase diagram and the Dirac point splits into 
pairs of Weyl points related by the unbroken symmetry. 
This emergent phase [7-21], dubbed as Weyl semimetal 
(WSM), is an appealing topological state of matter, with 
the Fermi surface being those Weyl points. 

In the simplest case when T symmetry is broken, a 
WSM at long wavelength only has one pair of Weyl 
points, which may be described by the Hamiltonian 

Pi = t[v x k x a x + VykyCJy + v z (k z - rh)cr z + b 0 \ • (1) 

Here Us are the Fermi velocities, cr’s are Pauli matri¬ 
ces, and r = =b denote the left- and right-handed Weyl 
fermions, which are required to come in pairs by the 
Nielsen-Ninomiya theorem [22]. rbz are the positions of 
the pair of Weyl points in the Brillouin zone (BZ), and 2bo 
is their energy splitting, which vanishes when V symme¬ 
try is not broken. Since all three cr’s are used up locally 
at each Weyl point, small perturbations may renormal¬ 
ize the parameters but cannot open a gap. Indeed, each 
Weyl point is protected by the Chern number (±1) of 
the valence band on a constant-energy surface enclosing 
it. Weyl points can only be annihilated in pairs of oppo¬ 
site chirality, when they are brought together (6, bo = 0) 
or when the translational symmetry is broken by strong 
interactions or by short-range scatterers. 

The topological properties of WSM can be best seen 
by considering a slice of the BZ normal to z. The Chern 
number of the slice changes from 0 to 1 and back to 0 as 
it crosses the two Weyl points successively. As a conse¬ 
quence, each nontrivial slice contributes an e 2 /h to the 
Hall conductivity [9-11] producing a xy = be 2 /irh in the 
bulk, and also contributes one edge state to the surface 
forming a surface Fermi arc connecting the two projected 
Weyl points. Recently, the Weyl points and surface arcs 
appear to be observed in optical experiments [23-25]. 


This progress may herald a flurry of exciting experiments 
on the appealing transport effects [26-28] predicted in 
WSM, e.g., the chiral magnetic effect [29-36] when bo be¬ 
comes nontrivial in the absence of V symmetry, and the 
axial magnetic effect [37, 38] when rb , viewed as a gauge 
field coupling oppositely to the left- and right-handed 
Weyl fermions, varies spatially. 

Here we discover a universal effect in lightly-doped 
WSMs, by examining the semiclassical dynamics of Weyl 
quasiparticles in the ballistic regime. The presence of 
Weyl points leads to substantial Berry curvatures. Un¬ 
der a longitudinal force field, the vector product of the 
force and the curvature results in a transverse velocity 
of the Weyl wave-packet, i.e., a velocity counterpart of 
the Lorentz force in the semiclassical description. There¬ 
fore, as the Weyl wave-packet propagates in the WSM, 
its center acquires a transverse shift upon transmission 
and reflection, perpendicular to the direction in which the 
Fermi energy or velocities change. The anomalous shifts 
are opposite for Weyl fermions with different chirality, 
and can be made imbalanced by breaking V symmetry. 
This effect has no counterpart in 2D and is analogous 
to the Imbert-Fedorov shift [39-44] observed for a circu¬ 
larly polarized light beam, a bosonic cousin of the Weyl 
wave-packet, due to variation of dielectric constant. 

We focus on the simple case of model (1) with bo = 0 
unless V symmetry breaking is noticed. Before the more 
general discussion based on the semiclassical equations of 
motion (EOM), we first derive our main results using the 
conservation of angular momentum J z , in the presence 
of rotational symmetry along z, which requires v x = v y 
and isotropic transverse perturbations. The total angular 
momentum of a Weyl wave-packet is given by 

J = r c x k c + ^n, (2) 

where (r c , k c ) is the phase-space center of the wave- 
packet, n is the unit vector (v x k x , v y k y , v z k z )/Ek, and 

Ek = rj^jv^k^PVykyPv^k^ are the energy dispersions of 

the electron and hole (rj = =t) bands. For simplicity, 
the subscript c will be dropped wherever appropriate. 
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The first term in J is the orbital angular momentum of 
the wave-packet, whereas the second term is the intrinsic 
1/2-pseudospin angular momentum of a Weyl fermion. 

We now consider the situation of a Weyl electron wave- 
packet impinging upon a potential step. As sketched in 
Fig. 1(a), consider an incident electron that has an energy 
Si > 0 and an incident angle 0\ = arctan(fc*/fc*) m th e x ~ 
z plane. The potential step satisfies V = 0 for z < 0 and 
V = Vo for 2 > 0 in our long-wavelength model, whereas 
it is smooth enough at the lattice scale. All energy scales 
are assumed to be small compared to the bandwidth such 
that the velocities change little across the step. 

For Vo < £i, the transmitted particle (T) is also a Weyl 
electron. The energy and momentum conservations re¬ 
quire that £i = Vo + St and k x == k x . The conservation 
of J z further determines the trajectory of the transmit¬ 
ted electron, regardless of the potential and scattering 
characteristics around z = 0. Remarkably, we find that 
in transmission the wave-packet center acquires a trans¬ 
lation perpendicular to the incident plane, given by 

x T _ r(nj-n^) _ tv z ( cot0 T cot#i^ ^ 
V ~ 2 U-Vo Si ) ’ ( J 

where = arctan(fcj/fcj) is the refraction angle. Since 

v z k J — \J {£i - Vo) 2 - v%k l x 2 , when V 0 > V c = £i~v x \k l x \, 
kj becomes imaginary and a total reflection occurs. 

For Vo > Si, the potential step forms a p-n junction 
and Klein tunneling of the Weyl electron may occur. 
The energy and momentum conservations require that 

k x = k l x and v z kj = —yj(Si — Vo) 2 — Evidently, 

when Si < Vo < 2 Si — V c , k J is imaginary and the total 
reflection occurs. Yet, when Vo > 2Si — V c , Klein tunnel¬ 
ing is possible and the minus sign in kj means a positive 
group velocity of the outgoing hole. Remarkably, the re¬ 
fraction index is negative, i.e, 6 l 6 T < 0. The transverse 
shift in this case takes the same form of Eq. (3). 

An anomalous transverse shift exists in reflection (R) 
as well. The conservation of J z directly leads to the trans¬ 
verse shift of the reflection plane by 

c R _ - n \) _ TV z cot 01 

V 2 k% Si ’ 1 ’ 

where k R = k l x and fc R = — k\. 

One observes that for an incident Weyl electron with 
r = 1, as shown in Fig. 1(b), and Sy T have the 
same sign that switches at 0i = 0. Sy T becomes peaked 
near total reflection. Approaching normal incidence, Sy T 
vanishes whereas fo/ R becomes divergent. The latter 
anomaly is resolved by the fact that the reflection prob¬ 
ability is exponentially small because of the pseudospin 
conservation. (Both reflection and transmission proba¬ 
bilities, determined by quantum mechanics, depend on 
Vo crucially.) As shown in Fig. 1(c), fo/ R is independent 
of Vo whereas Sy T strongly depends on Vq. Sy T is positive 
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FIG. 1. (a) Schematic figure showing the transverse shifts 

5y T,K in transmission and refection for an incident Weyl elec¬ 
tron wave-packet in the x-z plane scattered by a potential 
step, (b) and (c) Sy T,R plotted as functions of the inci¬ 
dent angle 0\ and the potential step height Vo, respectively. 
The shaded regions mark the range of total reflection, (d) 
Schematic figure showing the Fermi surfaces and momenta 
(k y — 0) of the reflected and/or transmitted wave-packets in 
normal, total reflection, and Klein tunneling regimes. The be¬ 
haviors in (a)-(c) are for left-handed Weyl electrons whereas 
the behaviors for right-handed ones (not shown) are opposite; 
the features in (d) do not distinguish chirality. We have cho¬ 
sen the parameter values as £\ — 0.01 eV, Vo = 0.014 eV, 
v z — 10 6 m/s, and v x , y = 3x 10 5 m/s in (b), and 0\ — 7r/4, 
£\ — 0.05 eV, v z — 10® m/s, and v XjV = 8 x 10 5 m/s in (c). 


for Vo < 0 and negative for 0 < Vo < V c . As Vo is further 
increased, the total reflection occurs until Vo = 2 S\ — V c , 
when the Klein tunneling begins to take place and then 
Sy T has the same sign as (Vo — 2Si). 

Physically, the shape of the Fermi surface changes 
across the potential step, as sketched in Fig. 1(d). Al¬ 
though k x is conserved, k z and hence the pseudospin 
angular momentum Tn z /2 changes during reflection or 
transmission. Therefore, the orbital motion must have 
a transverse shift to compensate this change, in order 





































3 


to guarantee the conservation of J z . As suggested by 
Eqs. (3) and (4), the transverse shifts are opposite for 
Weyl electrons with different chirality hence may be 
dubbed as a chirality Hall effect (CHE). 

While the symmetry argument based on J z conserva¬ 
tion can be applied to an interface that is sharp com¬ 
pared to the Fermi wavelength, the semiclassical descrip¬ 
tion of EOM is more general, regardless of any symmetry. 
Generalizing this formalism [45-48] to the lightly-doped 
WSMs, we find the following velocity and force EOM for 
the Weyl wave-packet center (r, fe): 

r = d k £ — ftkr • r — k x J7 , (5) 

k = —d r V — d r £ + Q r k • k . (6) 

We have included two possible spatially varying ingre¬ 
dients: the electrostatic potential V and the Fermi ve¬ 
locities Vi. The latter not only gives rise to the spatial 
dependence of the wave-packet energy £, but also leads to 
the phase-space Berry curvatures Q k r for the local eigen¬ 
states, where Q k . r . =d ki A rj -d rj A ki and A qi =i(u\d q .u) 
(q = r, k). In Eq. (5) 17 is the more familiar momentum- 
space Berry curvature with = e % i i Q kikj /2, and for the 
WSM described by model (1) 17 = —rv x v y v z k/(2£ 3 k ). 
The validity of the semiclassical description requires V 
and n’s to be slowly varying spatially compared to the 
Fermi wavelength. In this EOM we have considered cor¬ 
rections up to first-order in spatial gradients. 

In the presence of rotational symmetry along 5, the 
EOM (5) and (6) are simplified to k oc z and r = d k £ — 
QkzZ — k x 17, where Q kz = — rj(dk^>)(d z cos 0)/2 and 9 
(</>) is the polar (azimuth) angle of n. A straightforward 
calculation confirms that J z = 0. This demonstration 
shows that the various Berry curvature terms in the EOM 
are crucial for ensuring the conservation of J z , and thus 
for understanding the origin of CHE. 

When the rotational symmetry is broken, the argu¬ 
ment of J z conservation is not applicable. However, one 
can still use the EOM to obtain the transverse shifts for 
the same potential step. The potential step produces a 
force k — —d z Vz and the Weyl point leads to a substan¬ 
tial curvature 17. As seen in Eq. (5), these two effects 
result in an anomalous velocity — k x 17 of the wave- 
packet center perpendicular to z. This anomalous veloc¬ 
ity is analogous to the Lorentz force in the semiclassical 
picture. By integrating Eq. (5), we find an anomalous 
shift 





k x 17 dt 


TVxVy{r^_ - n\) 

2(vlkl + v%k$) 


(—k y , k x , 0), (7) 


where a = T for Vo <V C = £\ — ^J^ x k x + and a = R 

for Vo > V c . Apparently, Eq. (7) produces exactly the 
same results as Eqs. (3) and (4) in the case when v x = v y 
and Vo < £\. We note, however, two important issues for 
this semiclassical result of CHE. First, the semiclassical 


trajectory is unique, i.e., transmission if Vo < V c and 
reflection if Vo > V c . Second, Eq. (7) cannot be applied to 
the Klein tunneling case, because the place where £ k = 0 
requires a non-Abelian treatment [46], which is beyond 
the scope here. 

When V symmetry is broken, the two Weyl points with 
opposite chirality would split in energy, as denoted by rbo 
in model (1). In this case, all above results can still be 
applied by replacing £ a (a = I, R, T) by £ a — rb q. The 
imbalanced populations for the left and right chirality 
lead to a charge Hall effect. It is even possible to filter 
one chirality in transmission using the total reflection. 

A spatial variation of the Fermi velocities n’s can also 
lead to a CHE for Weyl wave-packets. Consider an in¬ 
terface around which n’s vary along z. In the presence of 
rotational symmetry, the conservation of J z dictates that 

5y a = r ^ 2fcI « cot 0 a - v l z cot 6>i) . (8) 

Interestingly, the total reflection would occur if e^) 2 > 
(' v l x ) 2 + (t^cot#i) 2 . In the absence of rotational sym¬ 
metry, one can still apply the EOM, assuming n’s to be 
slowly varying. For the outgoing electron, the anomalous 
shift in the x-y plane can be expressed as 

5r a = [ —4— [O kz d kz £ + (£2 x z)d z £] dt , (9) 

where whether a = R or T is definite and determined 
by the EOM, as in the case of Eq. (7). One observes 
that there is an additional transverse shift entirely due to 
apart from the shift due to the — k x J7 contribution 
that has been discussed above. Compared to momentum- 
space curvatures i7, the phase-space curvatures f2 rk are 
less well-known. Our result in Eq. (9) identifies a physical 
effect induced by 

The velocity-variation-induced CHE is similar to the 
Hall effect observed for circularly polarized light beams, 
i.e., Imbert-Fedorov shifts [39-44] due to the gradient 
of dielectric constant in isotropic media [41] in which 
Q k r = 0. Our result implies that Q kr can also be re¬ 
sponsible for a Hall effect of light beams in anisotropic 
media. Like photons, the Weyl fermions in WSM also 
have gapless and linear low-energy dispersions. However, 
their pseudospin is 1/2, the Lorentz invariance is broken 
at the lattice scale, their energies can be negative allow¬ 
ing Klein tunneling, and the separation of r = ± Weyl 
points in momentum is required to avoid mutual anni¬ 
hilation. In the reflection of close-to-normal incidence, 
the Imbert-Fedorov shift vanishes since the spin of light 
rn remains the same, whereas the CHE becomes anoma¬ 
lously large, as we have explained, since the pseudospin 
of Weyl electron rn/2 becomes opposite. As great ad¬ 
vantages due to their fermionic nature, gate voltages can 
easily modify the energetic landscape that Weyl fermions 
travel through; breaking V symmetry can easily create 
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FIG. 2. (a) Top view of an electron undergoing multiple (to¬ 
tal) reflections in a cylindric potential well of WSM. (b) Side 
view of the enhanced CHE in (a), (c) The chirality accumula¬ 
tion on the top and bottom surfaces in transmission through 
the green interface, which is detectable by the imbalanced 
absorbance of the left and right circularly polarized lights. 


imbalanced populations of the two Weyl flavors. Inter¬ 
estingly, for photons vacuum is a medium with the largest 
velocity, whereas for electrons vacuum can be viewed as a 
WSM with the strongest “intervalley scattering” leading 
to an infinite energy gap. However, as long as the surface 
Fermi arc appears, the CHE is anticipated to survive at 
a WSM-vacuum interface, with a total reflection. 

The CHE differs from the anomalous Hall effect (AHE) 
and the valley Hall effect (VHE). To view their dis¬ 
tinctions, we now reproduce the Hall conductivities of 
WSM using the EOM. An electric field produces a force 
k = —eE, leading to an anomalous velocity r = eE x n. 
Prom the current I = — e^ fc k, we derive that 

aij = eij *in J fr sgn( ™ A) • (10) 

For the model (1) with two Weyl points at rbz, a yz and 
g zx vanish due to the cancelation of Hall contributions 
from the two chirality r = ±, or from the opposite mo¬ 
menta ±k£. Yet, a xy survives because sgn (rv z rb) = 1. 
We can further obtain <j xy = be 2 /irh, assuming that the 
“gaps” rv z k z are “inverted” at 0 < rk z < b and taking 
into account the lattice regularization of the linear bands. 
For the same reasons the VHE, which requires an extra 
r-factor in the integrand of Eq. (10), simply vanishes. 

Evidently, in the above case the AHE arises from fi z 
(i.e., fi z leading to a xy = be 2 /nh and f2 x ( y ) producing 
a yz(zx) = 0): whereas the CHE can arise from any com¬ 
ponent of 17. Our results highlight that 17 is a two-form 
that has three components in the 3D case, instead of one 
component like the 2D case. In fact, the CHE has no 
counterpart in 2D, e.g., in graphene, because the inci¬ 
dent plane coincides with the atomic layer and a trans¬ 
verse shift normal to the layer is not well-defined. But an 
AHE (VHE) is possible in graphene when an energy gap 
opens at its Weyl points with T (V) symmetry breaking. 

Experimentally, the adiabatical variation of Fermi en¬ 
ergy and velocities may be controlled by local gates and 
strain effects, respectively. In a cylindric potential well, 
the incident Weyl wave-packets undergo multiple reflec¬ 
tions, leading to an enhanced CHE, as shown in Fig. 2(a) 


and Fig. 2(b). Because of the CHE in the multiple (total) 
reflections or in transmission, the left- and right- handed 
Weyl fermions are separated toward opposite surfaces, as 
sketched in Fig. 2(c). The surface chirality accumulation 
can be detected by the imbalanced absorbance of the left 
and right circularly polarized lights. When V symmetry 
is further broken, the imbalance between the two chiral¬ 
ity leads to a charge Hall effect with a voltage difference 
between the opposite surfaces. 

Finally, our predictions can be generalized to a sys¬ 
tem with multiple pairs of Weyl points [9, 12, 19, 24, 25]. 
Magnetic field effects can further be studied by including 
the Lorentz force eE x r and the orbital magnetic mo¬ 
ment m = e!7£ [46], and second order corrections can 
also be incorporated into the EOM [49] . Appealingly, the 
coupling E • B leads to the chiral anomaly [47], whereas 
the coupling —m-B shifts the Weyl wave-packet energy 
depending on its chirality. So far we have assumed that 
r6, the Weyl point locations, is insensitive to the spa¬ 
tial variation of V or u’s. If rb does change notably, 
V x b(r) would be a nontrivial axial gauge field, and we 
would have to take into account the corresponding axial 
magnetic effect. 

Note added .—After this work was finalized, a comple¬ 
mentary and independent study [50] appeared, with a 
similar CHE in reflection derived by a different approach. 
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